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13^ Abstract 

j-jTj The complex distance function Q , which plays a prominent role in the 

i— i definition of scalar (acoustic) wavelets, is found to determine a complex 

extension of the spherical coordinate system that is ideally suited for 
the construction of highly focused electromagnetic beams with helici- 
£T) ties conforming to the oblate spheroidal geometry of £■ This is used to 

t— I build a basis of electromagnetic wavelets F m = E m + iB m radiated or 

Q\ absorbed by the branch disk T> of Q. F m has integer angular momen- 

tum m around the axis of T> and definite spheroidal helicity. We use 
ON a regularization method to compute its singular charge-current den- 

sity and show that the total charge vanishes. Hence F m is due solely 
to electric and magnetic polarization currents. T> acts as a magnetic 
• • dipole antenna, and its axis as a coupled electric dipole antenna. We 

_ ^ propose this as an idealized electromagnetic model for quasars (with- 

^ out gravity, i.e., in flat spacetime), with V representing the accretion 

5_i disk and the vortex singularities along its axis representing the jets. In 

the regularized version, the accretion disk is represented by a solid, flat 
oblate spheroid and the jets by two solid, narrow semi-hyperboloids, 
as shown in Figure 2. 
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1 The complex helicity basis 

A rich geometry ensues when the Euclidean distance in R n is continued 
analytically [KO]. Here we specialize to M 3 and use the notation introduced 
in [HK9]. Define the complexification ( of the radial coordinate r by 

C = V( r - ia ) 2 = vV + z 2 = £ - in, r = (x, y, z) G R 3 , (1) 
a = ax, p = \J x 2 + y 2 , z = z — ia. 

The singularity of r at the origin expands to the branch circle of £, 

C = {r : ( = 0} = {(x, y, z) : x 2 + y 2 = a 2 , z = 0}, 

and continuation of £ along any simple closed loop linking with C gives — £• 
Hence £ is double- valued. To make it single- valued, it is necessary to close C 
with any membrane whose boundary is C, the simplest of which is the disk 

V = {(x, y, z) : x 2 + y 2 < a 2 , z = 0} = {r : £ = 0}. 

T> is a branch cut of and we choose the branch with £ > which reduces to 
+r when a — > 0. On crossing P, r] changes sign. As shown in the Appendix 
(48), (£,??) are related to the cylindrical coordinates (p, z) by 

az = £q, a 2 p 2 = (a 2 + £ 2 )(a 2 - r/ 2 ), 
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hence —a<r/<a. Far from T>, we have 

r 3> a £ ~ r and 77 ~ a • r = a cos 0. (2) 

The level sets of £ and 77 define an oblate spheroidal coordinate system, as 
explained in Figure 1 (see [HK9] for details). 




Figure 1: The real and imaginary parts of £ = £ — %r\ form an oblate spheroidal 
coordinate system with the z axis along a. Shown are cut-away views of an oblate 
spheroid £ with £ = e (e = 0.7a), a semi-hyperboloid H+ with 77 = 6 (b = 0.8a), and 
a semi-hyperboloid with 77 = —0' (6' = 0.5a). Also shown is the branch circle 
C = dT), which is spanned by the branch disk T>. C is the common focaJ circle of 
all the spheroids and hyperboloids. As e — > 0, £ shrinks to a double cover of T>. 
As b — > a, H± shrink to the positive and negative z axis. According to (2), the 
spheroids £ are deformations of the spheres r — constant and the semi-hyperboloids 
H± are deformations of the cones a cos 8 = constant, whose vertex r — expands 
to the circle {r : p = \J a 2 — r/ 2 , z = 0}. 

We now use £ to defined a complexified version •& of the polar coordinate 9: 

sint9 = -, cos$ = -. (3) 

Like i? is singular on C and discontinuous on the interior of T>. By applying 
the gradient operator with respect to r to the coordinates (£, </>), we obtain 
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vectors (£, 4>) extending the spherical basis (f , 6, 4>): 

% x 7 *•* r — ia pp + zz „ . „ 

C = VC = = = P sin # + z cos i? 

c c 

V-tf = — where i9 = p cos $ — z sin $ 



p C sin $ 

These vectors form a complex-orthogonal, oriented basis in C 3 , i.e., 



C X 1? = <j>, X <f> = C, x C = I?- 
In this basis, the gradient operator is 

v = c> + + t^^. 
C C sm ^ 

We now define a related basis which will be especially useful for studying 
electromagnetic radiation. Let 

/ n I N 1, 1 — COS l9 1, C — Z 

k = lntan(tf/2) = - In = - In > 



2 1 + costf 2 ( + z' 



so that 



Let 



d/t 



(4) 



d$ sin -i? 



1 

ip± = (f>± in, d± = -^j- = - (5^ T i<9 K ) 
and define the vectors 

„ , 4> : i 4>±i<d 

x± =V?P ± = ^±—- = ^ . 5 

p smif ( p 

Note that the factor d/t/di? transforms the singularity l/( of Vi? into the 
singularity l/p of V0, so that and 1? can be combined into ± This 
combination will be the key to helicity in our system, and (CiX+iX~) will 
be called the helicity basis. 
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Applying the chain rule in the coordinate system {^,tp^.,ij}—) gives the gra- 
dient operator in the helicity basis: 

V = Cd( + X+d+ + X-d- ■ 

Note that 

2 

C • X± = X± = 0, X+ ■ X- = -2 ( 6 ) 

C x x± = ±iX±, X+x-X- = —2, 

so x± are a P a i r °f null, transversal eigenvectors of the operator £x with 
eigenvalues ±i. We shall see that the eigenvalue i corresponds to advanced 
(absorbed) solutions and —i corresponds to retarded (emitted) solutions. 

In addition to the complex angular coordinates ip± , define the complex re- 
tarded and advanced time coordinates 

u = t — C, v = r + where t = t — is 

is complexified time [Kll] and we have chosen units in which the speed of 
light c = 1. The functions (u, v, ■0+, form a nuii coordinate system in 
spacetime related to the Newman-Penrose formalism. 

Note that (u,v) are singular on T> and ip± are singular on the z axis (5). 
Hence our coordinate system is singular on the set 

Q = {r:£ = 0}U{r : p = 0} = £> U U (7) 

where Z + and i?_ are the positive and negative z axis, respectively. We 
shall see that these coordinates are particularly well suited for representing 
electromagnetic fields radiated or absorbed by T>. The coordinate singular- 
ities will be natural for describing sources on Q, just the singularity of the 
spherical coordinate system (r, 9,(f>) at the origin is natural for describing 
point sources at r = 0. 



2 EM fields in the helicity basis 



We consider the class of real electromagnetic fields (E, B) represented in 
the complex coordinates (u,v,ip±) by 

E + iB = F = aC + p X+ + lX-, (8) 
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where 



2 2C 

is the complexified Coulomb field [N73] and a, f3, 7 are analytic functions of 
(u,v,ip±) when r ^ Q. The advantage of the representation (8) is that 

VxC = Vxx ± = 0, 

which will make it easy to compute VxF, and 

r i Q V • C = V • x± = 0, (9) 

which will make it easy to compute V • F outside of Q. 
The charge-current density (g,J) of F is defined 1 [K3] by 

g = V • F and — J = d t F + iV x F, (10) 

where the derivatives must be interpreted in a distributional sense and will 
be seen to be distributions (generalized functions) supported on Q. We shall 
compute these sources by regularizing the field F and then taking the limit 
as the regularization is removed. But first, let us find the conditions on 
(a, /3, 7) which make F a solution of the homogeneous Maxwell equations 
outside the singularity set Q. 

Note that 

dta = a u + ct v and d^a = a v — a u , (11) 

where the subscripts denote partial derivatives, and the same goes for f3, 7. 
Taking the divergence of (8) and using the properties of C and x± gives 



.a n ~ a u P- + 7+ 
C 2 + P 2 



2^p^ + 2 ^I±. (i 2 ) 



1 We use Lorentz-Heaviside units, where £0 = /io = 1, and set c = 1. The real parts of 
(10) give the inhomogeneous Maxwell equations with ( Re q, Re J) as the electric charge- 
current density, and imaginary parts state that ( Im g, Im J) is the magnetic charge- 
current density. To avoid magnetic monopoles, we will show that the total magnetic 
charge vanishes, so (Im g, Im J) is due solely to magnetic polarization. 



6 



Furthermore, 



d t F = 2 



C + (/?u + Pv)X+ + {lu+lv)X- 



, 2C 2a+ 2a- 
iVa x C = i{a + x+ + ot-X-) x ^2 = ~7^2~*+ ~ ~T2~X- 



x x+ = i{Pv ~ Pu)C x X+ + iP-X- ><X+ = (Ai ~ Pv)x+ + 



P 



27+ 



iV7 x x- = i(lv ~ 7u)C X X-+ il+X+ * X- = (lv ~ lu)X- ^C- 

P 



Collecting terms, we have 
-J = 2C 



a u + a v (3- - 7+ 



c 2 



r 



+ 2%, 



a 4 

C 2 



+ 2%. 



7v - -72 



c 2 



• (13) 



Requiring the sources (12) and (13) to vanish outside the singularity set Q, 
the homogeneous Maxwell equations are thus equivalent to 

(a v - a u ) sin 2 1? + /3_ + 7+ = 
(a u + a v ) sin 2 1? + /?_ - 7+ = 
a + - C 2 /? u = 



a- + C 7* = °- 



By simplifying and using 
1 



sin $ 



, cost? = tanh-^, ip = ^ + 1" — , C 2 = ^ 



cosh -0 ' 2^1 
this can be expressed in terms of of the coordinates (u, v, ip±) as 



7 + cosh 2 ip, = — /3_ cosh 2 ip, a + = ( 2 (3 U , a_ 



(14) 



These equations, which relate the longitudinal component a of the field to its 
transversal components /3, 7, express the homogeneous Maxwell equations in 
the coordinate system (u, t>,?/>±). 

Newman's holomorphic Coulomb field, representing the electromagnetic part 
of a spinning, charged (Kerr-Newman) black hole {i.e., ignoring gravity) 
[N73], is obtained by choosing 



1 



a 



, /? = 7 = 



C 2 



-V- 



,1 
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As proved in [K4], the charge-current density in this case is that of a disk 
T> spinning rigidly at the extreme relativistic angular velocity u = c/a, so 
that the boundary C = dT> spins at the speed of light. This kind of singular 
behavior is expected of black holes. 

The coherent electromagnetic wavelets constructed in [Kll], 

F = g'(r-C)X-, (15) 

are obtained by choosing 

a = /3 = 0, j = g'(u), 

which solves (14). Again, these are extreme solutions of Maxwell's equa- 
tions whose energy propagates at the speed of light everywhere outside the 
singularity set Q; see [Kllb, K12] for a detailed discussion. 

The function g' in (15) is the complex derivative of the analytic signal g of 
an arbitrary pulse function go(t), defined by 

1 [°° g (t')dt' 

This is a special case of the multidimensional analytic-signal transform [K3, 
Klla]. If go(t) is at all reasonable (e.g., square- integrable) , then g(t — is) is 
analytic outside the real axis and its real and imaginary parts are smeared 
version (to scale \s\) of go(t) and its Hilbert transform, respectively. The 
farther we get from the real time axis (by increasing |s|), the more smeared 
out and weaker the analytic signal g(t — is) becomes, and the same goes for 
its complex derivative g'(t — is). Since 

g'(i~ - = g'(t -£-i(s- 7/)), 

it follows that g'(r — C) is analytic for s > a (since then s — r] > s — a > 
for all r) and s < —a (since then s — r] < s + a < for all r). If s > a, 
g(r — Q peaks when rj = a (the positive z axis Z+), and if s < —a, it 
peaks when r\ = —a (the negative z axis ZJ). Hence g(r — Q is a beam- 
shaping pulse function, forming a beam along Z+ or Z- (depending on the 
choice of s) through its dependence on £. Such pulses are the basis for 
complex-source pulsed-beams in the engineering literature; see [HF89] and 
the detailed references given in [HK9]. 
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3 Electric-magnetic dipoles and quasar engines 



Equations (14) also have solutions representing non-axially symmetric gen- 
eralizations of (15). The simplest ones are 

F+ = (3(u, v, tP±)x+ F_ = 7 (u, v, ^±) X - (17) 

where (14) require that 

P u = /3_ = and 7^ = 7+ = 0, 

respectively. Thus F + must be advanced (absorbed) and independent of 
?/>-, and -F_ must be retarded (emitted) and independent of 

F + = /3(r + C, and F_ = 7 (r - C, 1>-)X- (18) 

where the dependence on <f> through ip± must be continuous and periodic. It 
is easy to verify directly that F± satisfy the homogeneous Maxwell equations 
outside Q, so their charge-current densities must be supported on Q. 

The fields (18) are null: 

F\ = (E± + iB±) 2 = 0, i.e., E\ - B\ = E± ■ B± = 0. (19) 

These are the conditions of pure radiation everywhere outside Q. As proved 
in [Kll, Kllb], this means that F± propagate without leaving any reactive 
(rest) energy behind; all their energy flows at the speed of light c, which is 
not true of generic EM fields, even in vacuum. Normally, electromagnetic 
energy flows at speeds v < c in the near zone due, among other things, to 
interference between field components propagating in different directions. 

Thus F± are extreme. To obtain 'normal' solutions whose energies flow at 
speeds less than c in the near zone, we must include a longitudinal compo- 
nent a(u,v,ip±) 7^ 0. This makes (14) more difficult to solve and will be 
pursued elsewhere. 

By (6), F± satisfies 

C x F ± = ±iF ± , (20) 

and we have seen that F± depends on time only through r ± We have 
interpreted this by saying that F + is absorbed and F_ is emitted. To 
confirm this, consider the limit of (6) as a — > 0, where £ — > f : 

rxF± = ±iF± ^ rxE± = ^B ± and f x B± = ±E±. (21) 
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This gives the Poynting vectors 

S± = E± x B± = =fE± x (f x E±) = =F#± f, (22) 

so the energy of F + flows inward to the origin and that of F- flows outward 
from the origin. This is indeed consistent with F+ being absorbed and F 
being emitted. 

Let us derive the spheroidal version of (22). It follows from (20) that the 
Poynting vector is given by 

T2S± = ±iFl xF± = F* ± x(Cx F±) = \F±\ 2 C - (C • F*±)F±, 

hence 

TS± = U±Z-±(Z-F*±)F± (23) 

where 

U ± = l -\F ± \ 2 = l -{El + B 2 ± ) = El 

is the energy density. The left side of (23) is real while the right side is 
complex. Since £ • F± = and C + C = 2V^, (23) can be rewritten as 

TS± = U±C - (V£ • F*±)F± = U±C - V£ • (E ± - iB±)(E± + iB±). 

Taking the real part of both sides gives 

TS± = E/±V£ - V£ • E±E± - V£ ■ B±B± = -Vf • T± (24) 

where T± is the Maxwell stress tensor of F± . This is the spheroidal version 
of (22), reducing to the latter as a — > since V£ —¥ f and • E±, V£ • B± 
both vanish. 

We now introduce the angular functions 

/ -i n \ m/2 

z rn = e -im^ ± = e ~im^ tan ±r»(0/ 2 ) = e~ ilmf> f + — ) (25) 

1 zt COS u 



< > ~ \ m/2 

<-"•'•' ( ± ' ]j with m c. 
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where m must be an integer for Z± to be single- valued. If ±m > 0, 
tan ±m ('[?/2) is a vortex factor suppressing the positive z axis and ampli- 
fying the negative z axis; if ±m < 0, it does the reverse. This can be seen 
by noting that z = £,r]/a (48), hence 

( ± z = — irj) ± (z — ia) = (£ ± z) =F i(a ± r/) = (a ± ry) 

a 

and 

C±z = 44> ?7 = =Fa «4> p = 0, =Fz > 0. 

Define the basic null helicity wavelets with angular momentum m about the 
z axis by 

F% = g'(T±()ZZx±- I (26) 



They form a Fourier basis for absorbed and emitted helicity fields analytic 
in ijj± and periodic in <p. 

To explain and justify the term helicity wavelets, start with the time-harmonic 
case. From (16) it follows that 

go (t) = e -^* g{t - is) = Sgns • H{ujs)e- i ^ t ~ ls \ (27) 

where H is the Heaviside function (see also [K3, page R298]). Hence 

g'(r±0 = -w Sgn(s±7 ? )^(o;(s±r ? ))e-^ ( ' ±f) e-^ ±,?) . 

At any fixed position r ^ Q, 

= Ae~ iujt X ± = - e" ia "(0 ± (28) 

where ^4 is a complex function of r, a, s, m and uj. If s > 0, then to > by 
(27) and (28) shows that the polarization of F™ spins clockwise about the 
vector r, as viewed from the origin, and that of -F™ spins counterclockwise. 
Since F- is retarded and F+ is advanced, both wavelets have negative 
helicity for ui > 0; that is, the polarization of F+ spins clockwise when 
viewed in the direction of propagation, toward the origin. If uj < 0, these 
spins are reversed. Thus we have shown: 

For time-harmonic signals, the wavelets F™ have helicity 2 — Sgnw. 



2 In [Kll] and [K4a], I used the opposite convention, so the helicity there was Sgnw. 
This can be arranged, for example, by defining r = t + is. 
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Furthermore, the above concept of helicity makes sense even when go is not 
time-harmonic. From (27) it follows that for a general signal, g{r) consists 
of the positive-frequency part of go if s > and the negative-frequency part 
of go if s < 0: 

1 f 00 

s > => g{r) = —J dco e-^-^oM = g^(r) 
S <0 9(r) = ~ f dcoe-^ s go(co) = -go(T), 
where <?o( w ) is the ordinary Fourier transform of go(t). 

• When s > a, F± is a superposition of negative-helicity fields. 

• When s < —a, F± is a superposition of positive-helicity fields. 

This idea extends to fields of indefinite helicity as follows. Define the helicity 
operator S on solutions as multiplication by — Sgnw in the Fourier domain. 
This can be expressed directly in the spacetime domain as —i times the 
temporal Hilbert transform: 

1 f'°° dt' 

SF = -j 7T-F(r, t') = F-(r, t) - F+(r, t), (29) 

where f denotes the principal-value integral and F^ denote the positive- 
and negative frequency parts of F. 

Remark 1. It may be thought unphysical to restrict the time-harmonic 
field (28) to just positive or just negative frequencies. This would be justified 
if the opposite component can be obtained by complex conjugation using a 
reality condition. 3 However, since F is complex, its positive- and negative- 
frequency components are independent. The reality conditions 

EUr)* = E_Ur), B w (r)* = B_ w (r) (30) 

apply only to the Fourier components of the real fields E,B. A general 
time-harmonic field has the complex form 

F(r, t) = F w (r)e"^* + F^(r)e iujt = E(r, t) + %B(r, t), u > 0, 

3 This is commonly done in the engineering literature, with the understanding that the 
true field is obtained by thing the real part. 
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where F w and F— u can be chosen independently. Then the real fields are 
E = Re F = ^ [F u + F*_Je~ iut + l -[F^ + F*Je iujt 
B = Im F = I[F W - Flje--* + I[F_ W - F*]e^, 

so 

as required. As we have seen, F w (r)e~ mi gives a real field of negative helicity 
and F -u{r)e lujt gives a real field of positive helicity. See [K4a]. 

Remark 2. Electromagnetic wavelets of definite helicity were already con- 
structed in 1994 [Klla]. However, those wavelets were globally sourceless, 
hence unrealizable. 



4 Regularized charge-current density of F± 

Let us write the general helicity fields (18) as 

F a (r,t) = g'(T + aC,^a)Xa, <r = ±l. (31) 

They satisfy the homogeneous Maxwell equations everywhere outside the 
singularity set Q: 

r <£Q V • F = and d t F a + iV x F a = 0, (32) 

hence their charge-current density (Q a ,J a ) must be a distribution supported 
on Q. We shall now compute this distribution by regularizing the field F a 
and then applying Maxwell's equations 'in reverse' to compute the sources. 
Let e > 0, < b < a, choose differentiable functions e(£) and h±(rj) with 

e«) = (°' « S ° M,) = ft "l" (33) 

and let 

R(£,v) = e{S)h+(T,)h-(-T,). (34) 

Then 

fl(0, ??) = R(£, ±a) = 0, hence r G Q => i2 = 0. 
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Furthermore, define the solid spheroid £ £ and semi-hyperboloids l~L h ± 

£ e = { r : < £ < e] H b ± = {r : b < ±77 < a} (35) 
and their union 

Q e ' b = £ £ UH b + UH b _. (36) 



Then 



r <£ Q £,b => R = 1. 



As e — > and b — > a, Q £ > — > Q. We shall regularize F a by multiplying it by 
R, thus eliminating the singularity on Q while leaving the field unchanged 
outside of Q £,b . We compute the charge-current density of the regularized 
field and then find the 'bare' charge-current density of F a as its limit. 

Thus define the regularized field 

F a = RF a . (37) 

Taking the divergence formally, we obtain the charge density of F a as 

~q g = V • F a = RV ■ F a + d^R V£ • F a + d v RVr) ■ F a 

= RV-F a + g'ir + a(, - id v )R, 

where we have used (50) and (51) from the Appendix. Assume that as 
t Qi R vanishes sufficiently rapidly that the first term is identically zero. 
(This simply requires some smoothness of R at £ = and rj = ±a.) Then 

Qa = ^^9 , (T + aC,ip a )d* c R, (38) 

where d£ = ^(d^ — id v ) is the partial derivative with respect to the complex 
conjugate £*. In terms of e nd h± we have 

Qa = Q% + Qt + Qa: ( 39 ) 

where 

g% = e'(Oh + (v)h-(-vM-Fa (40) 
gt = ±e(0ti ± (±r ] )V V -F a 
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and we have used the fact that 

h+(rj)^Q b<r]<a => h-{—rj) = 1 
h'-{—rj)^0 -a<r]<-b => h + (i]) = l. 

Similarly, we compute the current density of F a as 

-J a = d t F a + iV x F a 

= R(d t F a + iV x F a ) + idt:RV£ x F a + id v RVr] x F CT . 

Again, the first term vanishes and we have 

J a = -id{:RV{, xF„ - id v RVr) x F CT (41) 

= J a + J CT + J CT 

where 

Jl = -ie'(Oh+(v)h-(-vM x *V (42) 
Jt = Tie{Oh' ± (±r l )Vr ] xF (T . 



Since the supports of e'(£) and h'±(r)) are [0, e] and [6, a], respectively, it 
follows that (§0,J^.) are supported in £ e and {g^,J a ) are supported in 7^. 

The regularized currents can be expressed compactly by defining 



• V£xx CT io-f. iarj a 2 + £ 2 a 2 + £ 2 

''^7C = Z H \ ~2 2P + \ ~2 2 



_ i V?? xx„ o"7/ A <t£ / a 2 — t/ 2 „ a 2 — r] 2 2 

v ° = ~*^T7 = ~ z " TV ^T^ p+ V ^F 0, 

where we have used the identity ap = ^Jm^m v (48). According to (40) and 
(42), the currents are given by 4 

K = KqZ and Jt = viot- (44) 

Note that 

v e a -v e a = \ and ■ = 1, (45) 

4 Although v+ = > recall that ?? in and J^T is restricted to b < ±77 < a. 
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the significance of which will be investigated elsewhere. 5 

Equations (43) contain a wealth of information. For example, <t>-v% becomes 
infinite as r — > Z± (rj — > ±a), showing that £ £ has vortices there. Since 
is real, it relates the electric current density to the electric charge density 
and the magnetic current density to the magnetic charge density. On the 
other hand, the z and p components of v e a are imaginary, so they convert 
electric to magnetic sources and vice versa. 

The unregularized or 'bare' sources are obtained by taking the limits e — > 
and b — > a. By (33), 

e -)■ => e(£) -> H(£) and b —?■ a h±(r/) — )• H(a — rj), 

where H is the Heaviside step function. Therefore 

e'(O^S(0, h' + (ri)^-5(a- V ), -hL(~v) % + v) 

restrict the sources (g a ,J a ) to V and Z±, confirming that the bare sources 
are indeed distributions supported on these sets. Equations (40) give, after 
some simplification, 

^ + a z 

where we have dropped the superfluous factors H(^) and H(a =F v)- 
The unregularized versions of v e a and are 

e iarjp + a4> ± ^ 

« ff = and t»^r = T^z, 

where we have used the fact that \/ a 2 — rj 2 = p on V. 

5 It is tempting to interpret 5 % and 5^ as 'charge-flow velocities,' in which case (45) 
would state that the charges on £ E and "H^ now a t the speed of light c = 1. However, 
such an interpretation is valid only if the associated densities have a definite sign. For 
example, equal and opposite uniform charge distributions flowing in equal and opposite 
directions give a vanishing charge density but double the current. This suggests a 'two- 
fluid' interpretation whereby the sources are split into ones with positive and negative 
charges. (A further challenge is presented by the fact that v% is complex.) 
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5 The total charge vanishes: no monopoles 



As noted above, the regualrized sources g a ,Ja are complex and their real 
and imaginary parts are the electric and magnetic charge-current densities 
for F a . The total charge Q(t) must be constant due to charge conservation. 
But if this constant were nonzero, it would have to be reai to avoid magnetic 
monopoles, and this is impossible due to analyticity in r, unless Q(t) = 0. 
In this section we prove that Q(t) does indeed vanish identically. It suffices 
to prove this for the regularized basic helicity wavelets 



F™ = fl(£, r,)F™ = i?(£, V )g'(r + aC)e- im ^Xa- (46) 

The volume element in oblate spheroidal coordinates is given by [K3] 

f 2 + n 2 
d 3 r = ^ 1 d£dr?d(ft, 
a 



d6e- im ^ = / d6e- im ^ ±mK = 2tt<5° 



and 



since k = lntan(i?/2) is independent of eft. Equations (38) thus give the total 
charge of F a as 

/'oo pa 

Q™(t) = 2na5° m / d£ / dng'ir + <r£ - iarj) [d t R - id^R] . 

JO J -a 

Integrating the first term by parts in £ and the second term in rj, we obtain 

/a roo 
d V [Rg']™ =0 -2-Kia6 m J Q d£ [RgV v= - a 

poo pa 

-4trt< / d£ / d V Rd* c g'(T + aC). 

JO J -a 

The first boundary term vanishes because R(0, rj) = and 

£ — > oo g'(j + er£ — iar\) — > 0. 

The second boundary term vanishes because i?(£, ±a) = 0. Finally, the last 
term vanishes because g'(r + cr£) is analytic in £ = £ — ir/. This proves that 

Q?(t) ^ (47) 
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as claimed. Thus F a has neither magnetic nor electric monopoles. 

Since the total charge vanishes, the fields F a are due to a combination 
of electric and magnetic polarizations. Hence the support Q e,b (36) of the 
charge-current density acts as a combined electric-magnetic dipole antenna 

*~ TY). — 777, ~ ~ , 

for receiving F, or emitting F_ . If £ £ is flat and are narrow, meaning 
that 

e ^ , a — b 

- < 1 and < 1, 

a a 

then the electric dipole is supported mostly on U l~L b _ and the magnetic 
dipole is supported mostly on £ £ . 

We propose this system as an electromagnetic model for quasar engines,® 
with £ e representing the accretion disk and H ± the vortex jets; see Figure 2. 
It is known that quasars radiate light with a high degree of helicity (circular 
polarization) [BF2, E3]. Furthermore, since quasars are the most distant 
visible objects in the universe, they must radiate extremely powerful and 
highly collimated beams. The radiated wavelets F™ in (26) and their regu- 
larized versions F™ derive their power and collimation from a combination 
of the beam-shaping properties of g'{r — Q, their angular momentum, and 
the associated vortex factor [(£ + z)/(C — z)} m / 2 in (25). 

Remark. The above is a highly idealized and simplified model for quasars. 
Among other things, it ignores the fact that the accretion disk and the jets of 
real quasars consist of plasmas and the jets are coupled to the differentially 
rotating accretion disk by magnetohydro dynamic equations, or preferably 
a relativistic version thereof. Note that the shape of e(£) is related to the 
density of £ £ , and the shapes of h±(n) are related to the densities of If 
the plasma dynamics can be expressed in our spheroidal coordinate system, 
as suggested by Professor EnBlin (private communication) , then the coupling 
between the accretion disk and the jets could perhaps be represented by a 
relation between e(£) and h±(±rj). This would be an interesting subject for 
future study. 

6 We are obviously ignoring gravity by working in flat spacetime. Our model has roughly 
the same relation to quasars as Newman's holomorphic Coulomb field [N73, K4] has to 
spinning (Kerr-Newman) black holes. 
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Appendix 



Here we derive various expressions needed for the above computations. From 
the definition (1) and (3) we have 

( 2 = r 2 - a 2 - 2iaz = (f - ir/) 2 £ 2 - rf = r 2 - a 2 , az = £77 (48) 

2 2 2 2 2 2 2(f2 2\ 1 4 f 2 2 

a p = a r — a z = a (£ — 77 J + a — q rj = m^m v 

where 

2 1 a2 22 

= a + 4 , Tt?7 = a — ?? • 

1 . j_ ._ r — ia pp + zz „ . A -2 

£ = V£ — ? V77 = — - — = = p sin v + z cos v Q = 1 

"& = x C = zsint? — pcos$ •& = 1. 

Taking the real and imaginary parts of £ and simplifying gives 

V£=^±^ = ^p+^z, where n, = ^, p = |£| 2 
p p a P 

£a-r/r €n v ^ VP . , ™»? 

V77 = = z p, where n„ = — . 

pap p 

It follows that 

|2 



|V£| 2 = and |Vr?|' 



n- 



'/• 



To compute the regularized sources in Section 4, we need expressions for 
V£ " Xo-i ^ ' X CT an( i V£ x Xcr; x Xai where a = ±. It follows that 



V£xVr? = - — 0, (49) 
Z 1 



and since V£ -0 = 0, we have 



1? x C - Vf x C 

V£-Xa = itfVf ■ ~ = i<rVf ■ ~ — - = -«r0 • — 50 

P P P 

1 V£ x Vry a 
= — <70 • = a— . 

P P> 



Furthermore, 



C-X CT = (Ve-iVr ? )-x CT = Vryx CT = -^-. (51) 

A 1 
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Equations (50) and (51) will be used in Section 4 to compute the regularized 
charge densities. To compute the regularized currents, we need the following: 



vex 4> = ^ z -^p 

u a 



V£ x = V£ x (4> x C) = |V«e| 2 = 




V?? x = Vr] x (0 x C) = -i\Vr]\ 2 4> = 



in v cf), 



hence 




p + ia — cp 
P 



(52) 



[i ap 



Vr? x Xa 



— z 

/i ap 



p + cr — 0. 
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Figure 2: The charge-current densities (40), (42) of the regularized fields F a (37) 
are distributions supported on the union Q £ ' b (36) of a solid oblate spheroid £ £ 
and two solid semi-hyperboloids T-L b ± (35). According to (47), the total electric and 
magnetic charges vanish, hence Q £ ' b acts as an electric-magnetic dipole antenna. 
H± support the electric dipole and £ £ supports the magnetic dipole, and the flow 
of electric and magnetic currents within and between T-L b ± and £ £ is determined by 
the analytic signal g(r) (16). In the 'bare' (unregularized) limit (e — > 0, b — >• a), 
£ b approaches the branch disk V and T-L b ± approach Z±, so Q £ ' b — > Q (7). As 
explained in the text, Q £ ' b could be used to model quasar engines, with £ £ and T-L b ± 
representing the accretion disk and the jets, respectively. 
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